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ABSTRACT 

Although many methods of detecting extra-solar planets have been proposed and successful 
implementation of some of these methods enabled a rapidly increasing number of exoplanet 
detections, little has been discussed about the method of detecting satellites around exoplanets. 
In this paper, we test the feasibility of detecting satellites of exoplanets via microlensing. For this 
purpose, we investigate the effect of satellites in the magnification pattern near the region of the 
planet-induced perturbations by performing realistic simulations of Galactic bulge microlensing 
events. From this investigation, we find that although satellites can often cause alterations of 
magnification patterns, detecting satellite signals in lensing light curves will be very difficult 
because the signals are seriously smeared out by the severe finite source effect even for events 
involved with source stars with small angular radii. 

Subject headings: gravitational lensing - planets and satellites: general 



1. Introduction 

Various methods have been proposed to search 
for extrasolar planets (exoplanets). These meth- 
ods include the pulsar timing analysis, direct 
imaging, accurate measurement of astromctric 
displacements, radial velocity measurement, plan- 
etary transit, and gravitational microlensing [see 
the review of Perryman (2000)]. Since the first 
detection of an exoplanet around the pulsar PSR 
1257+12 (Wolszcan & Frail 1992), nearly 100 exo- 
planets have been identified (http:/ /exoplanets. org), 
mostly by the radial velocity method (Mayor & 
Queloz 1995). 

However, little has been discussed about the 
method of detecting satellites around exoplanets. 
This is mainly because it is thought to be too pre- 
mature to detect satellites given the difficulties of 
detecting exoplanets. Currently, the only promis- 



ing technique proposed to detect satellites is the 
transit method, where satellites are detected ei- 
ther by direct satellite transit or through perturba- 
tions in the transit timing of the satellite-hosting 
planet (Sartoretti & Schneider 1999). 

In this paper, we investigate the feasibility of 
detecting satellites of exoplanets via microlens- 
ing. Detection of a low-mass companion by us- 
ing microlensing is possible because the compan- 
ion can induce noticeable anomalies in the result- 
ing lensing light curves (Mao & Paczynski 1991; 
Gould & Loeb 1992). The microlensing method 
has an important advantage in detecting very low- 
mass companions over other methods because the 
strength of the companion's signal depends weakly 
on the companion/primary mass ratio although 
the duration of the signal becomes shorter with the 
decrease of the mass ratio. Then, if lensing events 
are monitored with a high enough frequency, it 
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may be possible to detect not only planets but 
also their satellites. Such a frequent lensing moni- 
toring program in space was recently proposed by 
Bennett (2000). 

The paper is organized as follows. In § 2, we de- 
scribe the microlensing basics of multiple-lens sys- 
tems. In § 3, we investigate the feasibility of satel- 
lite detections by carrying out realistic simulations 
of Galactic bulge microlensing events caused by an 
example lens system having a planet and a satel- 
lite. We summarize the results and conclude in § 
4. 

2. Multiple-lens Systems 

To describe the lensing behaviors of events 
caused by a lens system composed of a primary 
with a planet and subordinate satellites, it is re- 
quired to have the formalism of multiple-lens sys- 
tems. If a source located at ^(£,77) on the pro- 
jected plane of the sky is lensed by a lens system 
composed of iV-point masses, where the individ- 
ual components' masses and locations are m, and 
rL,i, the positions of the resulting images, n, are 
obtained by solving the lens equation, which is ex- 
pressed by 
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Here 6*e represents the angular Einstein ring ra- 
dius, which is related to the total mass of the lens 
system, m = m ii an d the distances to the 
lens, Dql, and the source, Dqs, by 
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The lensing process conserves the surface bright- 
ness of the source. Then the magnification of each 
image equals to the surface area ratio between the 
image and the unmagnificd source and mathemat- 
ically its value corresponds to the inverse of the 
Jacobian of the lens equation evaluated at the im- 
age position nj, i.e. 
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Although the magnifications of the individual im- 
ages cannot be measured due to the small sepa- 
rations between the images, one can measure the 



total magnification, i.e. A = Aj, where Ni 
is the total number of images. Note that to find 
the image positions and the magnification, it is 
required to invert the lens equation. 

For a single point-mass lens (N = 1), the lens 
equation can be easily inverted. Solving the equa- 
tion yields two solutions of image positions and 
the total magnification is expressed in a simple 
analytical form of 
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where u = (rs — »*l)/#e is the dimensionless lens- 
source separation vector normalized by #e- For a 
rectilinear lens-source transverse motion, the sep- 
aration vector is related to the single lensing pa- 
rameters by 
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where t E represents the time required for the 
source to transit 6*e (Einstein time scale), (3 is 
the closest lens-source separation in units of 6*e 
(impact parameter), t is the time of the maxi- 
mum magnification, and the unit vectors x and y 
are parallel with and normal to the direction of 
the relative lens-source transverse motion, respec- 
tively. The light curve of a single point-mass lens 
event is characterized by its smooth and symmet- 
ric shape (Paczyhski 1986). 

If the lens system has additional components 
(N > 2), the lens equation cannot be algebraically 
inverted. However, the lens equation can be ex- 
pressed polynomial in r\ and the positions 
of the individual images are obtained by numer- 
ically solving the polynomial (Witt 1990). If the 
lens system is composed of two lenses (e.g. the pri- 
mary and the plant), the lens equation is equiva- 
lent to a fifth-order polynomial in r\ and there 
exist three or five solutions of the image positions 
depending on the source location with respect to 
the lens components. The main new feature of the 
multiplc-lcns system is the caustics, which repre- 
sent the set of points in the source plane where 
the magnification of a point source becomes infin- 
ity, i.e. I det J| = 0. Hence a significant planet- 
induced deviation in the lensing light curve oc- 
curs when the source approaches the region around 
the caustic although the duration of the deviation 
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is short due to the small mass ratio between the 
planet and the primary, q p . The size of the caus- 
tic, and thus the probability of planet detections, 
also depends on the primary-planet separation and 
is maximized when the separation (normalized by 
Qe) is in the range of 0.6 <; d p <; 1.6 (lensing zone). 

As the number of lens components increases, 
solving the lens equation becomes nontrivial be- 
cause the order of the polynomial increases by 
N 2 + 1. One other method commonly used to 
obtain the magnification patterns of multiple- 
lens systems is the inverse ray-shooting technique 
(Schneider & Weiss 1986; Kayser, Refsdal & Sta- 
bell 1986; Wambsganss 1997). In this method, a 
large number of light rays are uniformly shot back- 
wards from the observer plane through the lens 
plane and then collected (binned) in the source 
plane. Then, the magnification pattern of the 
lens system is obtained by the ratio of the surface 
brightness (i.e., the number of rays per unit area) 
in the source plane to that in the observer plane. 
Once the magnification pattern is constructed, 
the light curve resulting from a particular source 
trajectory corresponds to the one-dimensional cut 
through the constructed magnification pattern. 
The advantage of using the ray-shooting method 
is that it allows one to study the lensing behav- 
ior regardless of the number of lens components. 
The disadvantage is that it requires a large com- 
putation time for the construction of the detailed 
magnification pattern. We attempted to investi- 
gate the lensing behaviors of multiple lens systems 
by solving the polynomial but we found that for 
the lens system of our interest for which the mass 
ratio between the least (satellite) and the most 
massive (primary) components is smaller than 
~ 10~ 5 , the numerical noise in the polynoimial 
coefficients due to the limited computer precision 
(~ 1CP 15 ) causes serious inaccuracy in solving the 
polynomial. We, therefore, use the ray-shooting 
method despite the requirement of large compu- 
tation time. 

3. Realistic Simulations 

To examine the feasibility of detecting satellites 
of exoplanets, we carry out realistic simulations 
of Galactic bulge microlcnsing events caused by a 
lens system having a planet and a satellite. The 
primary of the tested lens system is assumed to 



have a mass of 0.3 M & by adopting that of a late- 
type main-sequence star, which is believed to be 
the most common type of lenses for events de- 
tected towards the Galactic bulge (Alcock et al. 
2000). For the planet and the satellite, we test 
an Earth-mass planet and a Moon-mass satellite. 
Then, the mass ratios of the planet and the satel- 
lite with respect to the primary are q p = 10~ 5 
and q s = 1.2 x 10~ 7 , respectively. The planet is 
assumed to be separated by d p = 1.3 from the pri- 
mary. To investigate the dependence of the mag- 
nification pattern on the satellite position with re- 
spect to the planet, we test four cases of satel- 
lite locations with different combinations of the 
planet-satellite separation, d s (normalized by the 
Einstein ring radius of the planet 9e, p = v^^p^e)' 
and the orientation angle, a (measured from the 
primary-planet axis). 

For a low-mass companion, the source size 
can have a significant effect on the shape of 
the companion-induced anomalies in lensing light 
curves (Bennett & Rhie 1996). For the construc- 
tion of light curves, we, therefore, test three dif- 
ferent types of source stars with angular radii 
of p* = 8±/9 E = 1.5 x 10" 3 , 1.1 x 10- 3 , and 
0.85 x 10~ 3 . These values correspond to the source 
star radii of F0, GO, and K0 main-sequence stars 
located at Dos = 8.5 kpc for an event with a lens 
located at the half way point between the observer 
and the source, i.e. I?ol/-Dos = 0.5. We note that 
although a M-type star has a smaller source ra- 
dius (p+ = 0.6 x 10~ 3 for M0 main sequence), 
we do not consider the case because the star will 
be too faint to be observed (I > 25.6 even with- 
out extinction). For comparison, we also present 
the light curves corresponding to a source with 

= 3.0 x 10~ 5 , which is equivalent to the radius 
of a white dwarf. Although these events would be 
rare, there is at least some finite probability that a 
caustic crossing of a white dwarf would occur. In 
addition, although the brightness of a hot white 
dwarf would be faint (V ;> 25), this is not much 
fainter than a typical K dwarf. Therefore, instead 
of presenting light curves of events involved with 
a hypothetical source star having an arbitrary ra- 
dius, we present the light curves of events involved 
with a hot white dwarf. 

In Figure 1, we present the magnification pat- 
terns (gray-scale maps in the right panels) around 
the regions of deviations induced by the planet 
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Fig. 1. — Magnification patterns of a lens system composed of a primary, a planet, and a subordinate 
satellite (gray-scale maps in the right panels) along with the corresponding geometries of the lens systems 
(left panels). For both the magnification maps and the lens system geometries, the coordinates are centered 
at the center of mass of the lens system and all lengths are normalized by 6>e- The circle in the upper 
left panel and those in the other left panels represent the Einstein rings of the primary and the planet 
(with a radius # E;P = ^/%0^), respectively. The planet is separated by d p = 1.3 from the primary. The 
planet/primary and satellite/primary mass ratios are q p = 1CP 5 and q s = 1.2 x 10~ 7 , which correspond to 
a Earth-mass planet and a Moon-mass satellite around a 0.3M Q star, respectively. The labels in each of 
the left panel represent the planet-satellite separation, d s (normalized by #e, p ), and the orientation angle, a 
(measured from the primary-planet axis). The solid curve in the magnification map represents the caustics. 
For comparison, we present the magnification pattern unperturbed by the satellite in the upper right panel. 
The three white dots in the upper right panel represent the source sizes (corresponding to those of FO, GO, 
and KO main-sequence stars) which are used to construct the light curves in Fig. 2. The white lines in the 
right panels represent the source trajectories of the events, whose resulting light curves are presented in Fig. 
2. 
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Fig. 2. — Variation of the lensing light curve 
anomalies induced by an Earth-mass planet hav- 
ing a Moon-mass satellite. In each panel, we 
present three different light curves involved with 
source stars having different angular radii of = 
9+/9e = 0.85 x 10~ 3 (short-dashed curve), 1.1 x 
1CT 3 (long-dashed curve), and 1.5 x 10 -3 (dot- 
dashed curve), which correspond to the source 
radii of KO, GO, and FO main-sequence stars for a 
Galactic bulge event caused by a lens with a total 
mass m = 0.3M Q and located at Dol/Dos = 0.5. 
For comparison, we also present the light curves 
corresponding to the source of a white dwarf with 
p* = 3.0 x 10~ 5 , which are expected with (solid 
curve) and without the satellite (dotted curve). 
The source trajectories responsible for the light 
curves are marked on the magnification maps in 
Fig. 1, where the corresponding panels of the 
source trajectories and the light curves are marked 
by the same panel number. The time is presented 
in days assuming that the event has an Einstein 
timescale of = 30 days. The reference of the 
time is arbitrarily set. 



along with the corresponding geometries of the 
lens systems (left panels). In Figure 2, we also 
present the light curves of events resulting from 
the source trajectories marked in the correspond- 
ing panels in Fig. 1. From the simulations, we find 
the following results. 

1. If the satellite is located within the lensing 
zone of the planet, i.e. 0.6 <; d s <; 1.6 1 , the 
planet-induced caustic shape and the magni- 
fication pattern around the caustics are al- 
tered by the satellite due to the interference 
between the anomalies induced by the planet 
and the satellite [see panel (2) of Fig. 1]. 

2. If the satellite-planet separation is larger 
than the upper limit of the lensing zone 
of the planet, i.e. d s ;> 1.6, the interfer- 
ence becomes negligible. Then, the result- 
ing magnification pattern is well represented 
by the superposition of those of the two bi- 
nary systems where the planet-primary and 
the satellite-primary pairs act as indepen- 
dent lens systems (Bozza 1999; Han et al. 
2001; Han & Park 2002; Rattenbury ct al. 
2002). Although the satellite is located be- 
yond the planet's lensing zone, the satellite- 
planet separation is generally much smaller 
than the separation between the primary 
and the planet, i.e. <i s #E,p <C d p 8^. There- 
fore, the additional deviations induced by 
the satellite is located near the region of 
planet-induced deviations. 

3. If the satellite-planet separation is smaller 
than the lower limit of the lensing zone of the 
planet, i.e. d s <; 0.6, the planet-satellite pair 
behavior as if they are a single lens compo- 
nent with a mass equal to the combined one 
of the planet and the satellite. Since satel- 
lites are generally much less massive than 
their planets, i.e. q p + q s <~ q p , and thus 
the lensing behavior in this case is hardly 
affected by the presence of the satellite. 

4. The size (area) of the satellite-induced per- 
turbation region is comparable or smaller 

For a system of the Earth and the Moon located at 
Dql/^OS ~ 0-5, the Earth-Moon separation is about 3 
times of the Einstein ring radius of an Earth-mass planet 
in the foreground of the Galactic bulge. 



5 



than the size of source stars that can be 
monitorred from followup lcnsing observa- 
tions. As a result, the detailed structure in 
the lensing light curves is seriously smeared 
out due to the finite source effect. We 
find that the finite source effect is so severe 
that the satellite signals in the light curves 
of all tested events are completely washed 
out even for events involved with KO source 
stars. 

4. Summary and Conclusion 

We have tested the feasibility of detecting satel- 
lites by using microlcnsing. For this purpose, we 
have investigated the effect of satellites on the 
magnification pattern near the region of planet- 
induced perturbations by carrying out realistic 
simulations of Galactic bulge microlensing events. 
From this investigation, we find that although 
satellites can often affect the maginification pat- 
terns, detecting satellite signals in the lcnsing light 
curves will be very difficult because the signals are 
seriously smeared out by the severe finite source 
effect. 
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